library(tidyverse)
library(DeclareDesign)
library(rdss)
library(scales)

source("code/declarations/declaration_10a.R")

sims <- 2000
bootstrap_sims <- 2000

diagnosands <- 
  declare_diagnosands(
    var_estimate = pop.var(estimate),
    mean_var_hat = mean(std.error^2)
  )

diagnosis_10a <-
  diagnose_designs(
    declaration_10a,
    sims = sims,
    bootstrap_sims = bootstrap_sims,
    diagnosands = diagnosands
  )

write_rds(diagnosis_10a, file = "diagnosis_objects/diagnosis_10a.rds")

